Indications of Neutrino Oscillation in a 250 km Long-baseline Experiment 
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The K2K experiment observes indications of neutrino oscillation: a reduction of flux together 
with a distortion of the energy spectrum. Fifty-six beam neutrino events are observed in Super- 
Kamiokande (SK), 250 km from the neutrino production point, with an expectation of 80.1^5'4- 
Twenty-nine one ring /i-like events are used to reconstruct the neutrino energy spectrum, which is 
better matched to the expected spectrum with neutrino oscillation than without. The probability 
that the observed flux at SK is explained by statistical fluctuation without neutrino oscillation is 
less than 1%. 

PACS numbers: PACS numbers: 14.60.Pq, 13.15.+g, 23.40.Bw, 95.55.Vj 



Recent atmospheric and solar neutrino data indicate 
the existence of neutrino oscillation and therefore the 
existence of neutrino mass jl], |^, ||. The zenith angle 
distribution of atmospheric neutrinos observed by Super- 
Kamiokande (SK) shows a clear deficit of upward-going 
muon neutrinos, which is well explained by two-flavor u^- 
v T oscillations with Am 2 around 3 x 10 _3 eV 2 , and sin 2 20 
close to or equal to unity. 

The KEK to Kamioka long-baseline neutrino oscilla- 
tion experiment (K2K) [|| uses an accelerator-produced 
neutrino beam with a neutrino flight distance of 250 km 



to probe the same Am 2 region as that explored with 
atmospheric neutrinos. The neutrino beam is produced 
by a 12 GeV proton beam from the KEK proton syn- 
chrotron. After the protons hitting an aluminum target, 
the produced positively charged particles, mainly pions, 
are focused by a pair of pulsed magnetic horns Q . The 
neutrinos produced from the decays of these particles are 
98% pure muon neutrinos with a mean energy of 1.3 GeV. 

This analysis is based on data taken from June 1999 to 
July 2001, corresponding to 4.8 x 10 19 protons on target 
(POT). The pion momentum and angular distributions 
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downstream of the second horn are occasionally measured 
with a gas-Cherenkov detector (PIMON) |6| in order to 
verify the beam Monte Carlo (MC) simulation and to 
estimate the errors on the flux prediction at SK. The di- 
rection of the beam is monitored on a spill-by-spill basis 
by observing the profile of the muons from the pion de- 
cays with a set of ionization chambers and silicon pad 
detectors located just after the beam dump. The neu- 
trino beam itself is measured in a set of near neutrino 
detectors (ND) located 300 m from the proton target. 
The measurements made at the ND are used to verify 
the stability and the direction of the beam, and to de- 
termine the flux normalization and the energy spectrum 
before the neutrinos travel the 250 km to SK. The flux at 
SK is estimated from the flux of the ND by multiplying 
the Far/Near (F/N) ratio, the ratio of fluxes between the 
far detector (SK) and ND, to that of the ND. 

Since both a suppression in the number of events and 
a distortion of the spectrum are expected for neutrinos 
which travel a fixed path length in the presence of oscilla- 
tions, both the number of observed events and the spec- 
tral shape information at SK are compared with expecta- 
tion. All of the beam-induced neutrino events observed 
within the fiducial volume of SK are used to measure the 
overall suppression of flux. In order to study the spectral 
distortion, 1 ring ^i-like events (lRn) are selected to en- 
hance the fraction of charged-current (CC) quasi-elastic 
(QE) interactions {y ^ + n — > fi + p). Only the muon is 
visible in these reactions since the proton momentum is 
typically below Cherenkov threshold. The energy of the 
parent neutrino can be calculated by using the observed 
momentum of the muon, assuming QE interactions, and 
neglecting Fermi momentum: 

Eiec = rnNEj, - ml/2 

m N -£ M + P M cos0 M ' 1 ' 

where mjv, E^, m^, and # M are the nucleon mass, 
muon energy, the muon mass, the muon momentum and 
the scattering angle relative to the neutrino beam direc- 
tion, respectively. 

The ND is comprised of two detector systems: a 1 kilo- 
ton water Cherenkov detector (1KT) and a fine-grained 
detector (FGD) system. The flux normalization is mea- 
sured by the 1KT to estimate the expected number of 
events at SK. Since the 1KT has the same detector tech- 
nology as SK, most of systematic uncertainties on the 
measurement are canceled. The energy spectrum is mea- 
sured by analyzing the muon momentum and angular 
distributions in both detector systems. The 1KT has 
a high efficiency for reconstructing the momentum of 
muons below 1 GeV/c, and full Air coverage in solid angle. 
However, the f KT has little efficiency for reconstructing 
muons with momentum above 1.5 GeV/c since they exit 
the detector. The FGD, on the other hand, has high 
efficiency for measuring muons above 1 GeV/c, and the 



two complementary detectors are thus able to completely 
cover the relevant energy range. 

In the 1KT analysis, a cylindrical fiducial volume of 
25 tons oriented along the beam direction is used. Event 
selection criteria for the flux normalization are the same 
as those in reference Q . Events which deposit more than 
~ 100 MeV of energy are used for the measurement of the 
integrated flux. The measurement has a 5% systematic 
uncertainty, of which the largest contribution comes from 
the vertex reconstruction Q . For the spectrum measure- 
ment, further cuts are imposed in order to select lR/i 
events which stop in the detector. Among the events se- 
lected for the flux normalization measurement, 53% of 
the events have 1 ring. The events with a muon exiting 
from the detector are effectively eliminated by requir- 
ing the maximum charge of any PMT to be less than 
200 photo-electrons; 68% of the 1 ring events remain af- 
ter this requirement. The largest systematic uncertainty 
for the spectrum measurement is an uncertainty on the 
energy scale. The energy scale is understood to within 
!J %, which is confirmed with both cosmic-ray muons 
and beam-induced ir°s. 

The FGD is comprised of a scintillating fiber and water 
detector (SciFi) 0, a lead-glass calorimeter (LG), and a 
muon range detector (MRD) In the FGD analysis, 
events containing one or two tracks with vertex within 
the 5.9 ton fiducial volume of the SciFi are used. The 
track finding efficiency is 70% for a track passing through 
three layers of scintillating fiber and close to 100% for 
more than 5 layers S . Three layers is the minimum track 
length required in this analysis. Events which have at 
least one track passing into the MRD are chosen in order 
to select ^-induced CC interactions. The momentum of 
the track is measured by its range through the SciFi, LG, 
and MRD, with accuracy of 3%. 

If the proton produced in the QE interaction has a 
momentum greater than 600 MeV/c, its track may also 
be reconstructed. In the case where a second track is 
visible, the kinematic information is used to enhance the 
fraction of QE events in the sample. Assuming QE inter- 
action, the direction of the proton can be predicted from 
the muon momentum. The QE enhanced sample is se- 
lected by requiring that the direction of the second track 
agrees with the prediction to within 25 degrees. Events 
where the direction of the second track differs from the 
prediction by more than 30 degrees are put into a non- 
QE enhanced sample. In the QE enhanced sample, 62% 
of the events are estimated to be QE events from the 
MC simulation. In the non-QE enhanced sample, 82% 
of events are estimated to come from interactions other 
than QE. The SciFi events are divided into three event 
categories: 1-track, 2-track QE enhanced, and 2-track 
non-QE enhanced samples. 

The 2-dimcnsional distributions of the muon momen- 
tum versus angle with respect to the beam direction of 
four event categories (the 1KT event sample and the 
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TABLE I: The central values of the flux re- weighting parame- 
ters for the spectrum fit at ND ($nd) and the percentage size 
of the energy dependent systematic errors on $nd, F/N ratio, 
and £sk- The re-weighting parameters are given relative to 
the 1.0—1.5 GeV energy bin. 
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FIG. 1: (a) The muon momentum distribution of the 1KT 
\Rfi sample, (b) the angular distribution of the 1KT lRfi 
sample, (c) the muon momentum distribution of the SciFi QE 
enhanced sample, and (d) that of the SciFi non-QE enhanced 
sample. The crosses are data and the boxes are MC simulation 
with the best fit parameters. The hatched histogram shows 
the QE events estimated by MC simulation. 

three SciFi event samples) are used to measure the neu- 
trino spectrum. A x 2 -fitting method is used to com- 
pare these data against the MC expectation. The neu- 
trino spectrum is divided into 8 energy bins as defined 
in Table [jj During the fit, the flux in each energy bin 
is re- weighted relative to the values in the beam MC. 
These weights are normalized relative to 1.00 for the 
E v = 1.0—1.5 GeV bin. An overall normalization is in- 
cluded as a free parameter in the fit. The parameter 
R n qe is used to re-weight the ratio between the QE and 
non-QE cross-section relative to the MC simulation. The 
systematic uncertainties of the ND are incorporated into 
the fitting parameters. They are the energy scales, the 
track finding efficiencies, and the detector thresholds. In 
addition, the spectrum measurement by P1MON is used 
as a constraint on the re-weighting factors. 

The value of X 2 is 227.2/197 d.o.f. at the best fit point. 
All the parameters including the detector systematics are 
found to lie within their expected errors. The best fit val- 
ues of the flux re-weighting factors are shown in Table |. 
The muon momentum and angular distributions of li?/i 
events in the 1KT, and the muon momentum distribu- 
tions of the 2-track QE enhanced and non-QE enhanced 
events in SciFi are overlaid with the re-weighted MC in 
Figure [I]. The fit result agrees well with the data. The 
errors of the measurement are provided in the form of an 
error matrix. Correlations between the parameters of the 
fit are taken into account in the oscillation analysis using 
this matrix. The diagonal elements in the error matrix 
are shown in Table |[ 

The uncertainty due to neutrino interaction models is 
separately studied. In QE scattering, the axial vector 



mass in the dipole formula is set to a central value of 

1.1 GeV/c 2 and is varied by ±10%. The axial mass 
for single pion production is set to a central value of 

1.2 GeV/c 2 and is varied by ±20% g§. This affects 
both the q 2 dependence of the cross-section and the total 
cross-section. For coherent pion production, two differ- 
ent models are compared: one is the Rein and Sehgal 
model [jy], and the other is a model by Marteau Jl2| . 
For deep inelastic scattering, GRV94 |f| and the cor- 
rected structure function by Bodek and Yang Q are 
studied. For the oscillation analysis the Marteau model 
and Bodek and Yang structure functions are employed. 
Varying the choice of models causes the fitted value of 
Rnqe (= 0.93) to change by ~ 20%. In order to account 
for this, an additional systematic error of ±20% on R nqe 
is added to the error matrix. The choice of models and 
axial mass does not affect the spectrum measurement it- 
self beyond the size of the fitted errors. 

The F/N ratio from the beam simulation is used to 
extrapolate the measurements at the ND to those at SK. 
The errors including correlations above 1 GeV, where 
the PIMON is sensitive, are estimated based on the PI- 
MON measurements. The errors on the ratio for E v be- 
low 1 GeV are estimated based on the uncertainties in 
the hadron production models used in the K2K beam 
MC 0. The diagonal elements in the error matrix for 
the F/N ratio are summarized in Table Gj. 

The events in SK are selected using the timing informa- 
tion provided by the global positioning system. Events 
detected in SK must occur within an expected beam ar- 
rival time window of 1.5 /xsec. In addition, the detected 
events must have no activity in outer detector, and have 
an energy deposit greater than 30 MeV with a vertex re- 
constructed within the 22.5 kiloton fiducial volume [0. 
This sample of events is referred to as the fully con- 
tained (FC) sample. The efficiency of this selection is 
93% for CC interactions. Fifty-six events satisfy the cri- 
teria. With the timing cut, the expected number of at- 
mospheric neutrino background is approximately 10~ 3 
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events. 

The expected number of FC events at SK without os- 
cillation is estimated to be 80.1^g' 2 . The correlations 
between energy bins from the spectrum measurement at 
the ND and the F/N ratio are taken into account in the 
estimation of the systematic errors. The major contri- 
butions to the errors come from the uncertainties in the 
F/N ratio ( t.t'o%) an d the normalization (5.0%), dom- 
inated by uncertainties of the fiducial volumes due to 
vertex reconstruction both at the 1KT and SK. 

A two flavor neutrino oscillation analysis, with v p dis- 
appearance, is performed by the maximum-likelihood 
method. In the analysis, both the number of FC events 
and the energy spectrum shape for events are used. 
The likelihood is defined as C — C norm x C shape- The nor- 
malization term C norm {N O0S , N exp ) is the Poisson proba- 
bility to observe N b s events when the expected number 
of events is N exp (Am 2 , sin 2 29, /). The symbol / repre- 
sents a set of parameters constrained by the systematic 
errors. These parameters are described in detail later. 
The shape term, C shape = P(E t ; Am 2 , sin 2 29, f), 

is the product of the probability for each lJfyi event to 
be observed at E™ c — Ei, where P is the normalized 
El ec distribution estimated by MC simulation and Nm^ 
is the number of lRfi events. 

In the oscillation analysis, the whole data sample is 



used for C r , 



N, 



bs 



56. The spectrum shape in 



June 1999 was different from that for the rest of the run- 
ning period because the target radius and horn current 
were different. The estimation of energy correlations in 
the spectrum at the ND and in the far/near ratio has 
not been completed for this period. Thus, data taken in 
June 1999 are discarded for £ s hape- The discarded data 
correspond to 6.5% of total POT. The number of IR/j, 
events observed excluding the data of June 1999 is 29, 
and the corresponding number of lJfyt events expected 
from MC simulation in the case of no oscillation is 44. 

The parameters / consist of the re-weighted neutrino 
spectrum measured at the ND ($nd), the F/N ratio, 
the reconstruction efficiency (csk) of SK for lR/i events, 
the re- weighting factor for the QE/non-QE ratio R nqe , 
the SK energy scale and the overall normalization. The 
errors on the first 3 items depend on the energy and 
have correlations between each energy bin. The diagonal 
parts of their error matrices are summarized in Table | 
as described earlier. The error on the SK energy scale 
is 3% Two different approaches are taken for the 

treatment of systematic errors in the likelihood. The first 
is to treat the parameters / as fitting parameters with 
an additional constraint term in the likelihood (method 
1)0. The other approach is to average the £(/) sampled 
over many random trials weighted according to the prob- 
ability density distribution of the systematic parameters 
/ (method 2) @. 

The likelihood is calculated at each point in the Am 2 
and sin 2 29 space to search for the point where the like- 




FIG. 2: The reconstructed E v distribution for lRfi sample 
(from method 1). Points with error bars are data. Box his- 
togram is expected spectrum without oscillations, where the 
height of the box is the systematic error. The solid line is the 
best fit spectrum. These histograms are normalized by the 
number of events observed (29). In addition, the dashed line 
shows the expectation with no oscillations normalized to the 
expected number of events (44). 



lihood is maximized. The best fit point in the phys- 
ical region of oscillation parameter space is found to 
be at (sin 2 29, Am 2 )=(1.0, 2.8 x 10~ 3 cV 2 ) in method 
1 and at (1.0,2.7 x 10~ 3 eV 2 ) in method 2. If the 
whole space including the unphysical region is consid- 
ered the values are (1.03,2.8 x 10~ 3 eV 2 ) in method 1 
and (1.05, 2.7 x 10~ 3 eV 2 ) in method 2. The results from 
two methods are consistent with each other. At the best 
fit point in the physical region the total number of pre- 
dicted events is 54.2, which agrees with the observation 
of 56 within statistical error. The observed E]f c distribu- 
tion of the lRjj, sample is shown in Figure ^| together with 
the expected distributions for the best fit oscillation pa- 
rameters, and the expectation without oscillations. Con- 
sistency between the observed and best-fit E£ ec spectrum 
is checked by using Kolgomorov-Smirnov(KS) test. A KS 
probability of 79% is obtained. The best fit spectrum 
shape agrees with the observations. 

The probability that the observations are due to a 
statistical fluctuation instead of neutrino oscillation is 
estimated by computing the likelihood ratio of the no- 
oscillation case to the best fit point. The no-oscillation 
probabilities are calculated to be 0.7% and 0.4% for 
method 1 and 2 respectively. When only normalization 
(shape) information is used, the probabilities are 1.3% 
(16%) and 0.7% (14.3%) for the two methods. Allowed 
regions of oscillation parameters are evaluated by calcu- 
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FIG. 3: Allowed regions of oscillation parameters. Dashed, 
solid and dot-dashed lines are 68.4%, 90% and 99% C.L. con- 
tours, respectively. The best fit point is indicated by the star. 

lating the likelihood ratio of each point to the best fit 
point, and are drawn in Figure |3|. Both methods give 
essentially the same results. In order to be conserva- 
tive, the result from method 1 is shown in the figure as 
it gives a slightly larger allowed region at the 99% C.L. 
The 90% C.L. contour crosses the sin 2 29 — 1 axis at 1.5 
and 3.9 x 10~ 3 eV 2 for Am 2 . The oscillation parameters 
preferred by the total flux suppression and the energy 
distortions alone also agree well. Finally, the uncertain- 
ties of neutrino interactions are studied using the same 
procedure as the spectrum measurement at the ND. It is 
found that the effects of the interaction model difference 
on all the results are negligible due to the cancellation 
caused by using the same models in both the ND and 
SK. 

In conclusion, both the number of observed neutrino 
events and the observed energy spectrum at SK are con- 
sistent with neutrino oscillation. The probability that the 
measurements at SK are explained by statistical fluctua- 
tion is less than 1%. The measured oscillation parameters 
are consistent with the ones suggested by atmospheric 
neutrinos. At the time of this letter the K2K experiment 



has collected approximately one-half of its planned 10 20 
protons on target. 
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